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ABSTRACT 

We study the formation of non-linear structures in Warm Dark Matter (WDM) mod- 
els and in a Long-Lived Charged Massive Particle (CHAMP) model. CHAMPs with a 
decay lifetime of about 1 yr induce characteristic suppression in the matter power spec- 
trum at subgalactic scales through acoustic oscillations in the thermal background. We 
explore structure formation in such a model. We also study three WDM models, where 
the dark matter particles are produced through the following mechanisms: i) WDM 
particles are produced in the thermal background and then kinematically decoupled; 
ii) WDM particles are fermions produced by the decay of thermal heavy bosons; and 
hi) WDM particles are produced by the decay of non-relativistic heavy particles. We 
show that the linear matter power spectra for the three models are all characterised 
by the comoving Jeans scale at the matter-radiation equality. Furthermore, we can 
also describe the linear matter power spectrum for the Long-Lived CHAMP model in 
terms of a suitably defined characteristic cut-off scale fcctn similarly to the WDM mod- 
els. We perform large cosmological iV-body simulations to study the non-linear growth 
of structures in these four models. We compare the halo mass functions, the subhalo 
mass functions, and the radial distributions of subhalos in simulated Milky Way-size 
halos. For the characteristic cut-off scale k cut = 51/iMpc _1 , the subhalo abundance 
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(~ 10 9 Af sun ) is suppressed by a factor of ~ 10 compared with the standard ACDM 
model. We then study the models with A; cut — 51, 410, 820/iMpc -1 , and confirm that 
the halo and the subhalo abundances and the radial distributions of subhalos are indeed 
similar between the different WDM models and the Long-Lived CHAMP model. The 
result suggests that the cut-off scale k cut not only characterises the linear power spectra 
but also can be used to predict the non-linear clustering properties. The radial distribu- 
tion of subhalos in Milky Way-size halos is consistent with the observed distribution for 
^cut ~ 50 — 800/iMpc -1 ; such models resolve the so-called "missing satellite problem". 

Subject headings: cosmology: theory - early universe - dark matter 



1. Introduction 



The precise measurement of the cosmic microwave bac kground fCMB) aniso tropics established 



the standard A + Cold Dark Matter (ACDM) cosmology (jKomatsu et al.ll201ll ). Observations of 



the large-scale structure of the Universe, such as the galaxy power spectra from the Sloan Digital 
Sky Survey (SPSS) also co n firmed its success i n predicting the large scale structures of the Universe 
(e.g. iTegmark etaD faooj ) : iReid et alJ ()201oh : IPercival etafl fcoioh ). 



The validity of th e ACDM model on the galactic and the subgalactic scales has long been 
caught up in debate. iMoore et al.l (119991 ) argue that the number of dark matter s ubhalos is 
10 — 100 times larger than the number of satellites observed around the Milky Way (IKravtsov 



2010) . The so-called "missing satel li te problem" has been revisited in a somewha t quantitative 



cont ext ( Boylan-Kolchin et al.l 1201 ll ; lLovell et al.1 120121 ; iBovlan-Kolchin et al.ll2012l ). For exam- 
ple, iBovlan-Kolchin et al.1 (|201lh argue that, in the ACDM model, ~ 10 most massive subhalos 



in a galactic halo are too concentrated to be consistent with the kinematic data for the bright 
Mil ky Way satellites. Also, observations of the rotation velocities of galaxies using the 21 cm line 
by iPapastergis et al.l (|201ll ) show that the abundance of galaxies with observed velocity width 
w = 50kms~ 1 is ~ 8 times lower than predicted in the ACDM model. 



It is often suggested that WDM models resolve the apparent problems on subgalactic scales (]Bode et al 



2001 



). WDM particles have non- negligible velocity dispersions, which act as an effective "pressure" 
of the WDM fluid. Essentially, the subgalactic-scale density fluctuations are suppressed. The re- 
sultant matter power spectrum is quickly reduced around the cut-off scale that is determined by 
the velocity dispersion. Motivated by th e recent interest in this problem, seyeral authors study 



the s tructure formation in WDM models (IDunstan et al.l l201ll : ISchneider et al.l 12012 : iMenci et al 
20121 ). 



Constraints on WDM models can be obtained from astronomical observations. Observations of 



Lyman-a forests are often used for the purpose (]Viel et al.ll2005l : iBovarskv et al.ll2009i ). Absorption 
features in quasar spectra reflect the number density of neutral hydrogen, from which we can esti- 
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mate the matter power spectrum along the line of sight, even at large wavenum bers k ~ 10 h Mpc 



WDM models have also interesting impl ications for the cosmic reionization (jBarkana et all 12001 



Yoshida et al.ll2003l ; lGao Sz Theuns!l2007l ) . The formation of the first objects, and hence the produc- 
tion of ionizing photons, are delayed in WD M models. On the ot her hand, WDM models could help 
the completion of the cosmic reionization. lYue <fe Chenl (|2012l ) suggest that the reduced number 
of subhalos in WDM models makes the recombination of ionized hydrogens inefficient and results 
in earlier completion of the cosmic reionization. It is clearly important to study the clustering 
properties in WDM models in both linear and non-linear evolution regimes. 

There are also renewed interest in particle physics. Several candidates for WD M are suggested 



in pa rticle physics models beyond the St andard Model, such as light gravitinos (I Kawasaki et al 



1997 ). sterile neutrinos (see lKusenkd ( 20091 ) for a review and references) and super WIMPs (jCembranos et al 
20051 ) . It is important to notice that WDM particles can be produced via different mechanisms. 
Nevertheless, the above constraints from astronomical observations are focused on a single quantity, 
e.g., the mass of WDM particle in a specific model. It is unclear if such constraints can be applied 
to WDM models with different production mechanisms. Detailed comparisons of a wide class of 
models are clearly needed. 

In this paper, we also consider a Long-Lived CHAMP model. Throughout this paper, we 
assume that CHAMPs have an elementary charge, either positive or negative. CHAMPs are gen- 
erally realized in models beyond the Standard Model of particle physics. One such example is 
a slepton, a superpartner of leptons in super symmetric models. Sleptons as the lightest super- 
symmetry particles (LSPs) are stable when R-parity is conserved. The abundance of such stable 



CHAM Ps, however, is severely constrained by the searches in deep sea water (see iBeringer et al 



(J2012J) for a review and references). CHAMPs can also be unstable; a CHAMP decay into neu- 
tral dark matter and other decay products including at least one charged particle. For exam- 
ple, the stau can be the next lightest supersymmetric particle (NLSP) when the gravitino is 
the LSP (jBuchmuller et al.l 120061 ) . It is well-known that CHAMPs could affect the big bang nu- 



cleos y nthesis (BBN) reaction r ates and thus change the abundance of light elements (IPospeloy 
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with a lifetime about 1 yr can act effectively as WDM through acoustic oscillations in the thermal 
background. We study the effect of the oscillations on the matter power spectrum. 

We calculate the linear evolution of the matter density fluctuations for the three WDM models 
and the Long-Lived CHAMP model. We show that the comoving Jeans scale at the matter-radiation 
equality characterises the linear matter power spectra in the three WDM models well. We use the 
obtained linear matter power spectra as initial conditions of A-body simulations to follow the 
non-linear evolution of the matter distribution. We compare the halo mass functions, the subhalo 
mass functions, and the radial distributions of subhalos in Milky Way-size halos to discuss the 
clustering properties in the WDM models and in the Long-Lived CHAMP model. We show that 
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these statistics are similar when the cut-off scale is kept the same. We find that the WDM models 
and the Long-Lived CHAMP model with the characteristic cut-off scale A: cu t ~ 50 — 800 h Mpc -1 
resolve the so-called "missing satellite problem" . 

The rest of this paper is organised as follows. In Sec. [21 we summarize three WDM models and 
a Long-Lived CHAMP model we consider. Then, we introduce the common cut-off scale k cut which 
characterises the linear matter power spectra in these models. In Sec.[3j after describing the details 
of A-body simulations, we show simulation results and discuss their implications. Specifically, we 
mention the similarity of these models with the same cut-off and the possibility that CHAMPs 
behave like WD Ms and resolve the "missing satellite problem". Finally, in Sec.lU concluding 
remarks are given. 



Throughout this paper, we take the cosmological parameters that are given in lKomatsu et al 



(|201lh as the WMAP+BAO+H Mean; W0n b h 2 = 2.255, C dm^ 2 = 0.1126, ft A = 0.725, n s = 
0.968, r = 0.088 and A^(jfeo) = 2.430 x 10~ 9 , while we replace the energy density of CDM fi CDM /i 2 
by the energy density of WDM f2\yDM^ 2 f° r the WDM models and by the energy density of neutral 
dark matter produced by the CHAMP decay for the Long-Lived CHAMP model. 



2. WDM models and Long-Lived CHAMP model 

In this section, we summarize three WDM models and a Long-Lived CHAMP model we con- 
sider in this paper. We describe production mechanisms of WDM particles in each model and show 
the exact shapes of the velocity distribution. In the following subsections, we focus on three WDM 
models to specify our discussion, although our results can be applied to any WDM models with the 
same shape of the velocity distribution. Then, we introduce the Jeans scale at the matter-radiation 
equality. The matter power spectra in the three WDM models with the same Jeans scale at the 
matter-radiation equality are very similar. Their initial velocity distributions affect the damping 
tail of the matter power spectra. We also describe the evolution of the linear matter density fluctua- 
tions in a Long-Lived CHAMP model. The matter power spectrum is truncated around the horizon 
scale at the time when CHAMPs decay. Interestingly, the resulting power spectrum appears similar 
to those in WDM models. 



2.1. Thermal WDM 

In this type of models, fermionic WDM particles are produced in the thermal background. 
They are decoupled from the thermal background as the Universe expands and cools. At the 
time of the decoupling, their momentum obeys the thermal distribution, that is, the Fermi-Dirac 
distribution. We consider the generalized Fermi-Dirac distribution, 

f b) = eP /T W l + i • a) 
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Here and in the following, p denotes the comoving momentum of WDM particles, and Twdm is the 
effective temperatu re that character ises the comoving momentum of WDM particle s. In the case of 
the l ight gravitino (|Dine et al.l Il996l ) and the thermally produced sterile neutrino (jOlive fc Turner 



19820 . Twdm relates to the temperature of the left-handed neutrino T u through the conservation 

1 /3 

of the entropy, Twdm = ( tt^ ) T, where g& ec is the effective number of the massless degrees 

\ ydcc / 

of freedom at the decoupling from the thermal background. Note that f3 determines the overall 
normalization of the momentum di stribution and = 1 in t he case of the gravitino and the 
thermally produced sterile neutrino. iDodelson &: Widrowl ()1994l ) propose the sterile neutrino dark 
matter produced via active-sterile neutrino oscillations. In this case, the active neu trinos in th e 
thermal background turn into the sterile neutrino via the coherent forward scattering (|Clindll992l ). 
The resultant momentum distribution of the sterile neutrino is given by the generalized Fermi-Dirac 
distribution (see Eq. (p])) with Twdm — T, and j3 oc O^M where 6 m is the active-sterile mixing angle 
and M is the mass of the sterile neutrino. 



2.2. WDM produced by the thermal boson decay 

There are models in which the Majorana mass of the sterile neutrino arises from th e Yukawa 



coupling Y with a singlet boson ( Shaposhnikov Sz Tkachev 2006 ; Petraki Sz Kusenkoll2008 ). In these 



models, the singlet boson couples to the Standard Model Higgs boson through an extension of the 
Standard Models Higgs sector. The singlet Higgs boson has a vacuum expectation value (VEV) of 
the order of the electroweak scale when the electroweak symmetry breaks down. When the sterile 
neutrino is assumed to be WDM with a mass of an order of keV, the Yukawa coupling should 
be very small Y ~ O(10~ 8 ). This small Yukawa coupling makes the singlet boson decay to the 
two sterile neutrinos when the singlet boson is relativistic and is in equilibrium with the thermal 
background. Here, it should be noted that the sterile neutrino model is one specific example. 
In WDM models, where relativistic bosonic particles in equilibrium with the thermal background 
decay into fermionic WDM particles through the Yukawa interaction, WDM particles have the 
same resultant momentum distribution (see Eq. ([21 below). The r esultant momentum distribution 
is obtained by solving the Boltzmann equation (jBovanovskvl |2008| ) , 



f(p) 



(p/Twdm) 1 / 2 



55/2(p/Twdm) 



(2) 



where 



9v{x) 



x e" 



n=l 



n 



(3) 



Here, we have ignored the low momentum cut-off that ensures the Pauli blocking, while it does not 

1 /3 

change our results. The effective temperature is given by Twdm — T v with the effective 

number of massless degrees of freedom at the production of the sterile neutrino g pio ~ 100. The 
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normalization factor j3 is determined by the Yukawa coupling Y and the mass of the singlet Higgs 
boson M, (3 oc Y 2 M~ l . The velocity distribution have an enhancement /b oc p -1 / 2 at the low 
momentum p/?wdm ^ 1 ; since the sterile neutrinos with lower momenta are produced by the less 
boosted singlet boson, the decay rate of which is larger due to the absence of the time dilation. 
This enhancement indicates the "colder" (than the thermal WDM) property of the sterile neutrino 
dark matter produced by the decay of the singlet heavy boson. 



2.3. WDM produced by the non-relativistic particle decay 



In this type of models, a non-relativistic heavy particle decays into two particles, one or both 
of which become WDM. Supersymmetric theories realize this type of scenarios e.g. when the 
LSP is the gravitino and the NLSP is a neutralino. The relic abundance of the NLSP neutralino is 



determined at the tim e of chemical decoupling by the standard argument (jGondolo fc Gelmini 



Griest fe Seckel 



1991 



1991 



Eventually, the non-relativistic neutralinos decay into LSP gravitinos that 
become WDM. The particles produced by the decay of the moduli fields and of the inflaton fields are 
another candidates of this type of WDM (|Lin et alJl200ll ; iHisano et al.ll200ll : iKawasaki et al J 12009 ; 
Endo et al.l 12009 : iTakahashil 120081 ) . When we assume the heavy particle decay s in the radiation 



dominated era, the momentum distribu tion of the decay products is given by (IKaplinghatl 120051 ; 
Strigari et al-lboOTl : koyama et alJhoill ). 



f(p) 



(p/Twdm) 



exp(-p 2 / T WDM) 



(4) 



where Twdm is given by Twdm = PcmO,(td)/a(to) with the physical canter-of-mass momentum 
P cm , the scale factor a(t) at the decay time and at the present time to- We have defined as 
H(t = td) = 1/2t where H(t) is the Hubble parameter and r is the lifetime of the heavy particle. 



2.4. Jeans scale at the matter-radiation equality 



Now, we introduce two quantities to characterise the property of WDM. One is the present 

'WDM I 
Pcrit l*=*0' 



energy density of WDM, Qwdm = £ T m \t=t ■ Throughout this paper, we assume WDM particles 



account for all of the dark matter, letting Owdm^ 2 = 0.1126. Another important physical scale is 
the comoving Jeans scale at the matter radiation equality i eq , 



kj = a 



a 2 



(5) 

t t C Q 



with the gravitational constant G. Here, pm is the matter density and a 2 is the mean square of 
the velocity of the dark matter particles (see Eq. ([5]) below). Dark matter particles with fcj ~ 
100 — 1000 Mpc -1 are usually called WDM and expected to resolve the "missing satellite problem". 
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We note that, in the present paper, we do not consider whether or not a particular set of 
^ wdm and kj is in a viable region of the respective model. One such example is the gravitino 
WDM, a representative of the Thermal WDM model (see subsection l2.ip . This model has only two 
parameters, the effective number of the massless degrees of freedom at the decoupling g^ ec and the 
gravitino mass m 3 / 2 , to set SIvvdm an d kj. When we assume kj ~ SOMpc" 1 , these two parameters 
are determined as gdcc — 1000 and rfi 3 / 2 — IkeV. The effective number of the massless degrees of 
freedom at the decoupling of the gravitino is at most gdcc ~ 200 in the Minimal Super symmetric 
Standard Model (MSSM), a nd hence, another mechanism suc h as entropy production is needed to 
explain the gravitino WDM (|lbe et alJhoill ; llbe fc Satolboid ). 



2.5. Linear matter power spectra and Normalized velocity distribution 



We follow the evolution of the primordia l adiabatic fluctua tions for the three WDM models by 
modifying suitably the public software, CAMB (ILewis et al.ll2000l). We adopt the covariant multipol e 



perturbation approach for the massive neutrino (|Ma fc Bertschingerlll995l ; lLewis fc Challinorll2002l ) . 



We replace the Fermi-Dirac distribution of the massive neutrino by the momentum distributions of 
the WDM models discussed above. Our approach is valid when the WDM particles are kinematically 
decoupled at the cosmic time of interest. The Jeans scale of interest is around kj ~ O(100) Mpc" 1 . 
The primordial fluctuation of this wavenumber enters the horizon at T ~ O(10) keV. In a large class 
of WDM models, WDM particles are kinematically decoupled before the QCD phase transition, 
Tqcd ~ 100 MeV, and thus our calculation is valid. 

For comparison, we calculate the normalized velocity distributions for the WDM models and 
the linear matter power spectra extrapolated to the present time z = 0. The results are shown in 
Fig. Q] for kj = 51Mpc _1 , which correspond to m 3 /2 — 2keV for the thermally-produced gravitino 
WDM. Here, we have defined the dimensionless matter power spectra as, 



A(fc) = ^k 3 P(k) 

with the matter power spectra P(k). 

The velocity distribution g(v) is normalized as follows: 

i 

dvg(v) = 1 , 



(6) 



(7) 



1 There are variants of WDM models in which WDM particles are produced by the non-relativistic particle decay at 
late epochs. The matter power spectrum could be affected if the parent particles decay around the matter-radiation 
equality. In this case, the matter density during the radiation-dominated era is mostly contributed by the parent 
(cold) particles, rather than by the decay products, and then the density fluctuations of the cold, neutral parent 
particles can grow logarithmically. Note that the effect is more pronounced if the parent particles are charged (see 
subsection l2.6[) . 
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Fig. 1. — The normalized velocity distributions (left panel) and the dimensionless linear matter 
power spectra (right panel) for the standard CDM model and the three WDM models with fcj = 
51 Mpc -1 . 



dv v 2 g(v) 



a 



(8) 



with the variance (second moment) of the velocity a 2 . Note that the first equation is normalized 
with respect to the present energy density of WDM, Qwdm, whereas the second equation relates 
the velocity variance to the comoving Jeans scale at the matter radiation equality fej given by 
Eq. ©. One can then expect that the power spectra for the WDM models with the same ^wdm 
and kj are very similar, as seen in Fig.[TJ There, we see differences between the WDM models only 
in the damping tail of the power spectra at k > 10 h Mpc -1 . 



2.6. Long-Lived CHAMP and Cut-off scale 



Sigurdson &; Kamionkowskil (|2004l ) formulate the linearized evolution equations for fluctua- 
tions in a Long-Lived CHAMP model. They show that the subgalactic-scale matter density fluc- 
tuations are damped via a mechanism called "acoustic damping". The comoving horizon scale at 
wh ich CHAMP decays determines the cut-off sc ale of the matter power spectrum, which is defined 
by (jHisano et alJbood : iKohri fc TakahashJl201oh 



aH 



*=T"Ch ' 



(9) 



where H is the Hubble parameter and rch is the lifetime of CHAMP. Smaller-scale density fluc- 
tuations with k > kch enter the horizon before CHAMP decays and can not grow due to the 
acoustic oscillations of CHAMP in the thermal background. On the other hand, larger-scale den- 
sity fluctuations with k < k^^ grow logarithmically even after entering the horizon due to the 
gravitational instability as the density fluctuations of CDM. We assume CHAMPs decay in the 
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radiation dominated era. Then the comoving horizon scale at t = rch is evaluated as, 



k Ch = 2.2 Mpc" 1 x — 
yr 



1/2 ' <?Ch \V4 



5.363/ 



(10) 



where gch is the effective number of massless degrees of freedom when CHAMP decays. 



We need to consider three physical processes for the CHAMP model. First, we describe the neu- 
tralization of CHAMP. A positively charged particle may become neutral by forming a bound state 
with an electron e. Its binding energy is, however, almost the same as the hydrogen, E^ e ~ 13.6 eV. 
Hence, the positively charged particle keeps charged until its decay, since we assume CHAMP de- 
cays in the radiation dominated era. A negatively charged particle may become neutral by forming 
a bound state with a proton p. Its binding energy E^ p ~ 25keV is almost m4 He /m p ~ 2000 times 
larger than E^ e and hence is expected to make the negatively charged particle neutral at T ~ 1 keV. 
However, Helium 4 He is produced through BBN, with which a negatively charged particle may form 
a binding state. Its binding energy E h 4^ c ~ 337 keV is almost (Zi^/Zp) 2 x ?7T,4 Hc /?n p ~ 16 times 
larger than E^ p . It should be noted that even a nega tively charged particle bound with a proton is 
wrested by 4 He through a charge-exchange reaction ( Kamimura et al.ll2009i ). Therefore, when the 
yield of CHAMP Y"ch (Y = n/s with the number density n and the entropy density s) is smaller 
than the yield of the Helium i4 He , almost every negatively char ged particle forms a bind ing state 
with a helium nuclei, which has one positive elementary charge (IKohri fc Takahashill2010l ). 



Second, the decay products of CHAMP may lead to energy injection to the thermal background. 
The r esulting injection energy density is constrained from the photodissociation of BBN (jKawasaki et al 
200lh and CMB y- and /z- parameters (|Hu fc Siiklll993h . However, models with CHAMP with al- 
most the same mass with neutral dark matter are not severely constrained by BBN nor by CMB. We 
focus on such an "unconstrained" model. Note that the small mass splitting ensures the relatively 
long lifetime of CHAMP and the "coldness" of neutral dark matter. 



Finally, CHAMPs are tightly coupled with baryons before its decay. ISigurdson &; Kamionkowski 

(|2004 ) assume 6>baryon = #Ch where 9 is the divergence of the fluid velocity. This approximation is 
valid when the Coulomb scattering between baryons and CHAMPs is efficient, i.e., CHAMPs and 
baryons are tightly coupled. However, the constraints from the Catalyzed BBN essentially allow 
only heavy CHAMP with mch ^ 10 6 GeV for rch ^ 10 3 sec. It is unclear if the Coulomb scattering 
between baryons and such heavy CHAMPs is efficient. We have calculated the scattering efficiency 
and found that the tightly coupled approximation is indeed valid throu gh the epoch of interest 
for mch < 10 8 GeV. The details of the calc ulations are found elsewhere ( Kamadal let all ). In the 
following, we adopt the formulation given bv lSigurdson &: Kamionkowskil (|2004l ) using the tightly 
coupled approximation between baryons and CHAMPs. 



We modify CAMB (ILewis et al.ll2000l ) to follow the evol ution of density fluctuations i n the Long- 
Lived CHAMP model. The basic equations are given in ISigurdson &: Kamionkowskil (|2004l ). We 
obtain the power spectra for several rchS (see Eq. (fTOl) ). We find that the CHAMP matter power 
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Fig. 2. — We plot the dimensionless linear matter power spectrum in Long-Lived CHAMP model (t ~ 2.5yr ). 
We compare it with the same dimensionless linear matter power spectra in the WDM models as in Fig.Q] 
The oscillation around k ~ 9 ft, Mpc" 1 is the imprint of the "acoustic damping" . 

spectrum is very similar to the WDM models, as seen in Fig.[2l when fcch is set such that 

fccut = h ~ 45 k C h ■ (11) 

Hereafter, we use k cu t defined in the above as a characteristic parameter of the models we consider. 
The corresponding lifetime of CHAMP is r^h — 2.5 yr in the figure. The imprint of the CHAMP 
"acoustic damping" on the linear matter power spectra is clearly seen. One can naively guess that 
structures in the Long-Lived CHAMP model would be similar to those in the WDM models. It is 
important to study the non-linear growth of the matter distributions in the Long-Lived CHAMP 
model. We use large cosmological A-body simulations to this end. 



3. Numerical simulations 



Our simulation code is the parallel Tree-Particle Mesh code, GADGET-2 (|Springelll2005l ). We 
use N = 512 3 particles in a comoving volume of L = 10/t Mpc on a side. The mass of a 
simulation particle is 5.67 x 10 5 h~ 1 M sun a nd the gravitationa l softening length is 1 h^ 1 kpc. We 
run a friends-of- friends (FoF) group finder ([Davis et al.1 Il985l ) to locate groups of galaxies. We 
also identify substru ctures (subhalos) in each FoF group using SUB-FIND algorithm developed by 
Springel et al.l ((20011). We do not assig n any thermal vel ocity to simulation particles because it can 
lead to formation of spurious objects ([Cohn et al,ll2008l ). We start our simulation from relatively 
low redshift z = 19, at which the thermal motion of WDM is redshifted and negligible. It should 
be noted that the heavy, neutral dark matter produced by the CHAMP decay is assumed to have 
negligible thermal velocities. 



In Fig. [31 we plot the projected matter distribution in the CDM model (left panel), in the 



Fig. 3. — The projected matter distribution in the CDM model (left panel), in the Thermal WDM model 
(middle panel) and in the Long-Lived CHAMP model (right panel). For the Thermal WDM model and for 
the Long-Lived CHAMP model, we take the same cut-off scale k cut — 51 Mpc" 1 as in Fig.\2[ One side of the 
plotted region is L = 10 /i" 1 Mpc. Brighter regions denote higher matter densities. 



Thermal WDM model (see subsection !2.ip (middle panel) and in the Long-Lived CHAMP model 
(right panel). For the Thermal WDM model and for the Long-Lived CHAMP model, we set the 
same cut-off scale k cut = 51 Mpc -1 as in Fig.[2j One side of the plotted region is 10 /i" 1 Mpc. 
Regions with high matter densities appear bright in the plot. We see that many small objects, i.e., 
halos and subhalos, have formed in the CDM model. Contrastingly, in the Thermal WDM model 
and in the Long-Lived CHAMP model, the matter distribution is much smoother and appears more 
filamentary. The abundance of small objects is much reduced. Overall, the matter distributions in 
the Thermal WDM model and in the Long-Lived CHAMP model look similar. Note that numerous 
small objects along the filaments in the Thermal WDM model and in the Long-Lived CHAMP 
model could be numerical artifacts; t his is a long-standing problem of hot /warm dar k matter 



simulations due to discreteness effects (jWang <fe White! 120071 ; iPolisenskv <fe Ricottil l201ll ) . Earlier 



studies propose a simple formula for the critical halo mass, 

M c = 10.1 x p M d mcSLn /c~ c 2 ak . (12) 

below which the abundance of halos is unreliable. Here, d mea n — L/N 1 ^ 3 is the mean comoving 
distance between simulation particles and fc pe ak is the wavenumber at the maximum of the A(k). 
We will discuss this point further in the following section. 

We compare the halo mass functions in the models we consider in Fig.^J The fiducial cut-off 
scale is k cut = 51 Mpc" 1 (left panel) as in Fig.[2j but we also show the results for fc cut = 410 Mpc" 1 
(right panel). The latter corresponds to m. 3 / 2 ~ 9.5 keV for the thermally-produced gravitino 
WDM (see subsection l2.ip and TQh — 0.04 yr for the Long-Lived CHAMP model. The halo mass 
M corresponding to k CVLt is given by 

M = *epf*y aax 10 , h .r M>im x (!iM 3 . (13) 

3 \k cut J V k cut J 
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Fig. 4. — The halo mass functions in the CDM model, in the three WDM models and in the Long-Lived 
CHAMP model with k cut = 51 MpcT 1 (left panel) and with k cut = 410 Mpc -1 (right panel). The upturn at 
the halo mass M ~ 2 x 10 9 h" 1 M sun (vertical line) in the left panel may be owing to the artificial objects 
due to the discreteness effects. 



distance from the center 


0-50 


50 - 100 


100 - 150 


150 - 200 


number of satellites 


19.7 


15.62 


8.08 


14.16 



Table 1: The number of ob served satellites in each 50 kpc from the center of our Milky Way. According to 
Polisenskv & Ricotti \201 we count observed satellites known before the SDSS as one and those found by 
the SDSS as 3.54 due to the limited sky coverage of SDSS. 



It is important to examine if the halo abundance is compromised by the above-mentioned 
numerical artifacts. For k cut = 51 Mpc -1 , we see upturns in the mass functions at M ~ 2 x 
10 9 h~ l M sun . This is owing to peculiar discreten ess effects in hot /warm dark matter simulations 
(|Wang fc White! 120071 : rPolisenskv Ricottil boill ). The critical halo mass (see Eq. (fl2"]) 1 is M c ~ 
2 x 10 9 /i" 1 M sun for our simulation parameters L = 10/i _1 Mpc, N = 512 3 and A: pea k — 3 Mpc" 1 
(see Fig. [2]). The estimated mass limit is indeed consistent with the upturn seen in the left panel 
of Fig.H 

Therefore, we conservatively restrict our discussion to halos with masses M > M c = 2 x 
10 9 h- 1 M SUQ for k cut = 51 Mpc" 1 . The number of halos at M ~ M c in the three WDM models 
and in the Long-Lived CHAMP model is ~ 10 times smaller than that in the CDM model. Note 
the similarity of the halo abundances in the three WDM models and in the Long-Lived CHAMP 
model, as naively expected from the similarity in the linear matter power spectra. 

In order to see if the WDM models and the Long-Lived CHAMP model resolve the "missing 
satellite problem", we select Milky Way-size halos with masses of 0.5 x 10 12 /i" 1 M sun < Afi ia i D < 
1.5 x 10 12 h -1 Msun i n our simulations. Although the halo mass of Milky Way itself is in debate 



(e.g. lXue et al.l (|2008l ) and references therein), we take a relatively lower mass among suggested 
values. For a (slightly) small value of Mh a i , the relative mass ratio M sate ii itc /Mh a i becomes larger. 
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Fig. 5. — The cumulative subhalo mass functions averaged over Milky Way-size halos in the CDM model, 
in the three WDM models and in the Long-Lived CHAMP model with fc cu t = 51Mpc _1 (left panel) and 
with fc cut = 410 Mpc -1 (right panel). The shaded region corresponds to the mass of nonlinear objects at 
which the upturn (numerical artifacts) occurs in halo mass function (see the vertical line in the left panel of 
Fig.\§. 



Then the apparent discrepancy between the number of observed satellites and the simulated subhalo 
abundance at a given mass scale becomes smaller. We thus choose the small Milky Way mass as a 
"conservative" one. 

We compare the cumulative subhalo mass functions averaged over the Milky Way-size halos 
in our five models in Fig.0 Note that we see again an upturn for fc cut = 51Mpc _1 around 
M su b/Mhost — 2 x 10~ 3 . Above the mass scale, where we can measure the mass function robustly, 
the subhalo abundance is suppressed by a factor of ~ 10 in the models with k cut = 51 Mpc -1 (left 
panel) compared with the CDM model. For the models with k cut = 410 Mpc -1 (right panel), the 
subhalo abundance is suppressed at most by a factor of ~ 2. 

Let us now examine the radial distribution of the subhalos in our simulated Milky Way-size 
halos. We adopt the abundance of observed satellites in our Milky Way listed in Tabled! We account 
for the sky coverage of SDSS as follows. We count the number of the observed satellites in each 
50 kpc bin such that each satellite known before the Sloan Digital Sky S urvey (SDSS) is weighted a s 
one whereas each satellite discovered by the SDSS is weighted as 3.54 (|Polisenskv k Ricottil I20I1I ) . 



In Fig.[U we compare the averaged radial distributions of the subhalos in the Milky Way-size 
halos. The left panel shows the radial distribution in the CDM model and the right panel shows 
the radial distribution in the Thermal WDM model with A; cu t — 810 Mpc" 1 , which corresponds 
to — 16keV for the thermally-produced gravitino WDM. We include all the subhalos with 

M > 2 x 10 7 h~ l M sim . We also show variation of the radial distribution by thin solid lines. It 
should be noted that we set the gravitational softening length to be 1 h^ 1 kpc. The subhalo count 
in the innermost bin could have been affected by the spatial resolution. We can see the CDM model 
predicts a larger number of subhalos by a factor of 2 — 10 in each radial bin than observed. Note 
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Fig. 6. — The radial distributions of subhalos in Milky Way-size halos in the CDM model (right panel) and 
in the Thermal WDM models (left panel) with k C ut — 810 Mpc -1 . We divide the distance from the center of 
host halo in 50 kpc bins. Each thin line corresponds to the radial distribution in each Milky Way-size halo. 
Thick lines represent the average over the Milky Way-size halos. For comparison, we also plot the radial 
distribution of the observed satellites listed in Tabic [I] 



also that the number of satellites roughly scales with the host halo mass. Among the five Milky 
Way-size halos we selected, which have masses of 0.5 x 10 12 h~ l M sun < M halo < 1.5 x 10 12 /i -1 M sun , 
the total number of subhalos differs by a factor of ~ 3. By comparing the two panels in Fig.[6j we 
find that the radial distribution of subhalos in the Thermal WDM model with fe cut — 810 Mpc -1 
is similar to the one in CDM model and hence, the Thermal WDM model with k cu t — 810 Mpc -1 
does not seem to resolve the "missing satellite problem". 

We plot the subhalo radial distribution in the Thermal WDM model with k cu t — 51 Mpc -1 in 
Fig. For this model, where the suppression of the subgalactic-scale structure is most significant, 
a sizeable fraction of subhalos have masses smaller than M c (see Eq. (|12|) ). Thus the number count 
in the radial distribution is likely unreliable. Note however that, even without discarding the small 
mass subhalos (M su b < M c ), the subhalo abundance is slightly smaller than the observed satellites. 

We plot the radial distributions of subhalos with masses M > M c ~ 10 8 h' 1 M sun in the 
Thermal WDM model with fc cut — 410 Mpc" 1 in Fig.[HJ For this model, discarding small subhalos 
(M su b < M c ~ 10 8 -Ms Un ) reduces the subhalo abundance by a factor of ~ 2. Overall, the radial 
distribution of subhalos (after discarding) appears to reproduce the observed distribution. 

Let us now examine closely the similarity of the three WDM models and the Long-Lived 
CHAMP model when the characteristic cut-off scale &; cu t is kept the same. We compare the averaged 
radial distributions of the subhalos in these models for k cu t — 410 Mpc -1 (again after discarding 
small subhalos) and for A; cut — 810 Mpc -1 in Fig.EO From Fig. HI Fig. [5] and Fig.EjJ we conclude 
that the similar cut-off scale in the linear matter power spectra yields also similar halo and subhalo 
abundances and radial distributions. 
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Fig. 7. — The radial distribution of subhalos in Milky Way-size halo in the Thermal WDM model with 
fccut = 51 Mpc~ . This subhalos may include the artificial small objects due to the discreteness effects. 

4. Summary 

In this paper, we study the formation of non-linear objects in three WDM models and in a Long- 
Lived CHAMP model. We calculate the time evolution of the matter density fluctuations in the 
linear evolution regime by suitably modifying the public software CAMB. By using the obtained linear 
matter power spectra as initial conditions, we also perform large cosmological iV-body simulations. 
The results are summarized as follows. 

First, the comoving Jeans scale at the matter-radiation equality characterizes the linear matter 
power spectra of WDM models well. In the three WDM models motivated by particle physics, 
WDM particles are produced in different ways, but the linear matter power spectra with the same 
Jeans scale are very similar except for some difference at the damping tail at large k. We also 
consider a Long-Lived CHAMP model which has been suggested to yield a cut-off of the matter 
power spectrum through the "acoustic damping" . The cut-off scale of the matter power spectrum 
in the Long-Lived CHAMP model is determined by the comoving horizon scale when CHAMP 
decays. We empirically find the correspondence of the cut-off scales k cut (see Eq. (jlip ) between the 
three WDM models and the Long-Lived CHAMP model. 

By performing large cosmological iV-body simulations, we compare the abundances of nonlinear 
halos and subhalos and the radial distributions of the subhalos in Milky Way-size halos. The 
three WDM models and the Long-Lived CHAMP model produce very similar halo and subhalo 
mass functions and radial distributions if k cu t is kept the same. Therefore, we conclude that A; cut 
determines the clustering property of WDM and Long-Lived CHAMP well in both linear and 
non-linear growth of the matter density. 

One might naively guess that our simulation results for small non- linear objects may be com- 
promised by numerical effects. However, our conclusions are drawn after discarding small objects 
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Fig. 8. — The same plot as Fig.\2\in the Thermal WDM model with k cut = 410 Mpc" 1 . After discarding 
small subhalos (M su b < M c ~ 10 8 h^ 1 M sun ), the subhalo abundance is reduced by a factor of ~ 2 (right 
panel). 



that are likely numerical artifacts. We also compare the subhalo radial distributions in Milky Way- 
size halos with that of the observed satellites. We find that the WDM models and the Long-Lived 
CHAMP model are broadly consistent with the observation when they have k cu t ~ 50 — 800 Mpc" 1 . 
This cut-off scale corresponds to m 3 / 2 ~ 2 — 16keV for the thermally- produced gravitino WDM 
and rch ~ 0.01 — 2.5 yr for the Long-Lived CHAMP. Because there is sig nificant variation of th e 
subhalo abundance among host halos with different masses, as reported by llshiyama et al.l (|2009i ). 
it would be important to use a large sample of halos in order to address the validity of the models 
in a statistically complete manner. 

Our results have a further implication for particles physics. We clarified how to put constraints 
on a few parameters of particles physics models which provide a WDM candidate. By calculating 
and comparing two quantities, the relic density of dark matter and the comoving Jeans scale at 
the matter-radiation equality, one can apply the reported constraints in a specific particle physics 
model (e.g. the mass of thermal l y-pro duced gravitino WDM or sterile neutrino WDM produced 
through the lDodelson fc Widrowi (119941 ) mechanism) to virtually any model parameters of interest. 



Finally, we note that the so-called astrop hysical feedback pr ocesses are also thought to affect 
the abundance of luminous satellite galaxies (jBrooks et al.l |2012| ) . Unfortunately, whether or not 
and how the baryonic processes changes our simple understanding based on the cut-off scale k cut ~ 
50 — 800 Mpc" 1 is unclear due to the complexity of the baryonic processes. Regarding "missing 
satellite problem", both WDM and Long-Lived CHAMP models and the baryon feedback may 
explain the small number of the observed luminous satellites. The degeneracy can be resolved, in 
principl e, by the direct prob es of non-luminous small objects, by, for instance, future submillilensing 
surveys (jHisano et al.1 120061 ) . 
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Fig. 9. — Comparison of the averaged radial distributions in the CDM model, in the three WDM models 
and in the Long-Lived CHAMP model for k cut ~ 410 MpcT 1 (left panel) and for k cut ~ 820 Mpc -1 (right 
panel). For k cut ~ 410 Mpc^ 1 (left panel), we have discarded small subhalos (M su b < M c ~ 10 8 h^ 1 M SVLn ) 
as in Fig.\S[ 
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